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Abstract The aim of this paper is to construct and analyze solutions to a 
^ , class of Hamilton-Jacobi-Bellman equations with range bounds on the optimal 

response variable. Using the Riccati transformation we derive and analyze a 
fully nonlinear parabolic partial differential equation for the optimal response 
function. We construct monotone traveling wave solutions and identify para- 
metric regions for which the traveling wave solution has a positive or negative 
wave speed. 
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■ 1 Introduction 

>: 

•/^ . The purpose of this paper is to analyze special solutions to a fully nonlinear 

partial differential equation which can be derived from the Hamilton-Jacobi- 
Bellman (HJB) equation for the value function arising in a class of optimal 
allocation problems. In many practical stochastic dynamic optimization prob- 
lems, the goal is to maximize the expected value of the terminal utility. More 
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precisely, let us suppose that X — Xf,t G [0,T], is a stochastic process satis- 
fying a stochastic differential equation (SDE): 

xt = ^,tix^)t + a!{xf)w„ (1) 

where /ij and af > arc the drift and volatility of Ito's stochastic process 
([T]). Here Wt {t > 0) denotes the standard Wiener process. The goal is to 
find an optimal response strategy {6} = {0t\t € [0,r]} belonging to a set A 
of admissible strategies and yielding the maximal expected utility from the 
terminal value X^, i.e., 

ma.xE\u(Xi)\X^ ^ x] . (2) 

{0}£A L V J /I u J 

In this paper we consider the case when the optimal response strategy 6 is 
restricted by the unity from above, i.e., A = {{0} \ 0t < 1, < t < T}. The 
function u represents the terminal utility function. Throughout the paper we 
shall assume that u is a strictly increasing and concave function, i.e., u'{x) > 
and u"{x) < for all x £R. 

It follows from the theory of stochastic dynamic programming (see e.g. [2]) 
that problem © can be solved by introducing the so-called value function 

V{x,t) := sup E[u{X^) I Xf = x]. (3) 
{9}eA 

Using the Bellman optimality principle combined with the tower law of con- 
ditioned expectations it can be shown that the value function satisfies the 
so-called Hamilton- Jacobi-Bellman (HJB) equation 

(4) 

for aU X e R and t g [0,r) (see e.g. PDHIISIE] ) . 

The main goal of this paper is to construct monotone traveling wave so- 
lutions to the HJB equation (|4]) subject to the constraint A = {{9} \ 6t < 
1, < t < T} for the optimal decision policy {6}. Depending on the models 
considered, we show the existence of traveling wave solutions with positive as 
well as negative wave speeds. 

This paper is organized as follows. In the next section we investigate a 
simple HJB equation with drift and volatility functions linearly depending 
on the optimal decision parameter 0. Using a Riccati-like transformation we 
transform the HJB equation, originally stated for the value function V into a 
fully nonlinear parabolic PDE for the reciprocal value of the optimal response 
function 9. We construct a traveling wave solution with a decreasing wave 
profile, and we extend the results to the case when the underlying processes 
is governed by a SDE with a drift quadratically depending on the parameter 
9. In section 3 we investigate a more general HJB equation with a volatility 
function depending nonlinearly on the optimal decision policy parameter 9. 
We again identify a range of model parameters for which a traveling wave 
solution exists and has a monotonically increasing profile. 
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2 Construction of a traveling wave solution to the HJB equation 
with a positive wave speed 

In this section, we focus our attention to traveling wave solutions to the HJB 
equation. First, we shall examine a simplified model where the drift fif{Xt) and 
volatility (jf{Xt) are linear functions in 9. Then we generalize the results to 
the HJB equation (j?]) with an underlying stochastic process satisfying a SDE 
(III) with a drift function quadratically depending on the control parameter 6. 



2.1 A simple HJB equation 

In what follows, we shall analyze a simplified model capturing essential 
features of a more complex model with general drift and volatility functions fi 
and (T. We assume is a Brownian motion with drift /i^ — utd and volatility 
a" = 0>O, i.e., 

xf = w6't + 6'^yt, 

where w > is a positive parameter. We restrict our strategy by 1 from above. 
Then the corresponding Hamilton-Jacobi-Bellman equation ([3]) for the value 
function V{x,t) reads as follows: 

— {x,t)+s^^^,^[-0' — {x,t)+ioe—{x,t)j^O, V{x,T)=u{x). (5) 

Suppose, for a moment, that ([5]) has a classical solution V such that dxV{x, t) > 
0, and d^V{x,t) < for all a; S M and t e [0,r]. For justification of such an 
assumption we refer the reader to Proposition [31 Let us denote by 9*{x,t) the 
optimal response strategy at {x,t). It maximizes the function 

1 d'^V dV 
m3 9^ -9^ — {x,t)+uj9—{x,t) G M 

subject to the constraint 9 < 1. li the optimal response satisfies 9*{x,t) < 1 
then we have 

n*f ,N d^V{x,t) 
Placing (jni) back into ([S]) we obtain 



dt 2 dlV 



= for < t < T, V(T, x) = u{x). (7) 



Following [T] and [3] (see also [S]), we introduce the Riccati-like transfor- 
mation 

, , ldlVix,t) 

tp(x,t) := ITTTT^. 8 

Recall that, in the context of optimal portfolio allocation problems, the func- 
tion (p is related to the so-called Arrow-Pratt coefficient of the absolute risk 
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aversion (c.f. [TU], [12 )■ Performing straightforward calculations, the evolution 
equation for Lp becomes 



dip 1 d r 1 dip\ 
'dt ^ 2'dx\'^'dx) 



= 0. (9) 



In view of the optimal response function 9* ~ l/ip^ equation ^ is fulfilled 
by (f in the region {(x, i), '^{x, t) > 1}. On the other hand, if ip{x, t) < 1, then 
the maximum in ^ is attained &t 9* = 1. Inserting 9 = \ into ^ we end up 
with an equation for V{x,t) of the form: 

dV 1 d'^V dV 

m^2^^-^-'^ nx,T) = .(.). 

In terms of the transformed function if, it can be further reduced to 

It is easy to see that equation (|10l) is satisfied in the region {{x, t), ip{x, t) < 1}. 

Combining equations ([9]) and (jlO[) allows us to rewrite them in a compact 
form ^ 

^ + + ^B{^) =0, a: G R, i G (0, T), (11) 

where 

i^, for^<l, r -^(1-^)2, for^<l, 

-2^, forv5>l, 0, for^>l. 

(12) 

Notice that A{if) and B{(f) are increasing and continuous functions. 
2.1.1 Construction of a traveling wave solution with a positive wave speed 

In this subsection we shall construct a traveling wave solution to pT|) of the 
form 

Lp[x,t) = v{x + c{T -t)), a: G M, t G [0,r], 

with the wave speed c G M where v = v{^) is a function defined on R. 

Inserting the aforementioned ansatz on the solution into ([TT|) we conclude 
the existence of a constant Kq G M such that the function v — v(^) fulfills the 
identity 

- cv{0 + ^AA{vm + B{v{i)) = ifo, (13) 
for all ^ G M. Let us introduce an auxiliary function z = z(^) as follows: 

2(0 = ^(t'(0)- 

Then <y9(a;, t) = v{x + c(T — t)) is a traveling wave solution to (jlip if and only 
if the function z is a solution to the ODE: 

z'iO = FiziO), (14) 
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Fig. 1 A graph of the function F{z) in the case c > and Kq + c < 0. The roots of 
the function F{z) (left) and those of the function F{z) = F{z) — z (right). 



where F{z) ^ Kq + cA-^{z) - B{A-~^(z)). Clearly, 

^,,(Ko + 2cz + ^il~2z)\ forz<i, 

Notice that the function F is continuous for z < 1. Its graph, for the case 
when c > and Kq + c < 0, is depicted in Fig. [1] In this case the function F 
has exactly two roots z^ such that F{z^) = and < z+ < 1/2 < z~ < 1, 
where 

1 c 1 



.- = 1 + ^, z-^-----Vcy^^-2ic + KQ)/.. (16) 

We have F'{z^) > and F'{z+) < and F{z) < for z+ < z < z^ . Therefore, 
up to a shift in the argument ^, there exists a unique solution z = z(^) to (jl4[) 
connecting the steady states z^ such that 

lim z{^) = z^, z+<z{^)<z~, z'(O<0, CeM. 

The traveling wave profile v ^ v{^) is given by = A-^{z{C)). It satisfies: 

lim = v+<v{^)<v~, < 0, C e M, 

with = A~^{z^). Clearly < v+ < 1 < . On the other hand, we 
can prescribe the limiting values and calculate the corresponding wave 
speed c > and constant Kq. Indeed, it is straightforward to verify that, for 
< w+ < 1 < , the wave speed c and Kq are given by formulae: 

c=-^ f, Kq = -cv-. 17 

Since is a smooth nonlinear function we obtain that the solution z 
is a smooth function in the ^ variable. However, as A is just smooth 
we obtain that the traveling wave profile f is a smooth function in the ^ 
variable only. As a consequence, we have that the solution ip(x,t) = v(x + 
c{T — t)) is smooth and it is a weak solution to (jlip in the usual sense. 



6 



Naoyuki ISHIMURA, Daniel SEVCOVIC 



2.2 A HJB equation for a quadratic drift function 

In this section we shall assume the underlying stochastic process satisfying a 
SDE ([T]) with a drift function /x^ quadratically depending on the parameter 6, 
i.e. := uiO — ^6^. The volatility is again assumed to be linear in 9, := Q. 
The Hamilton- Jacobi-Bellman equation Q for the value function V has the 
form 



'dt 



(18) 

In this case, the optimal response strategy d*(x,t) is given by 
provided that 9*{x,t) < 1. The function V solves the nonlinear PDE: 

'm~Y dW~d,.v ^^ forO<i<T, y(r,x) = x gm. (20) 

Applying the Riccati-like transformation 

l d^,V{x,t)-d,V{x,t) 1 ( dlV{x,t) \ 

ip{x, t) := = - 1 , 21) 

uj d^V{x,t) uj \dxV[x,t) J 

it is easy to verify that the equation for (p reads as follows: 

^ + A±^0 (22) 

dt 2 9x V^2 Q^) dx2g> ' ^ ' 

and it is satisfied by tp in the region {{x,t), 'fiix,t) > 1}. If ip{x,t) < 1, then 
the maximum in is attained at 6** = 1. In this case, the equation for the 
solution V{x, t) and its Riccati transformation (p are as follows: 

dV Id'^V 1 dV 

and 

provided that (p[x,t) < 1. Hence we can rewrite the equation for (p as follows: 

^ + ^A(^) + - 0, X- G R, t G (0, T), (24) 

where A{ip) = A{tp) and B{(p) = B{if) + A{(p) {A and B are defined in (jl2p ). 

Next, following the analysis from section [^.1.11 we can construct a traveling 
wave solution to ([M)) of the form 

^{x,t) ^i{x + c{T -t)), a; G M, < G [0,T], 
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with the wave speed c S M and the profile v = v{^). Since A = A and B = B+A 
the transformed wave profile z(^) = A{v{£,)) should satisfy the ODE: 

-z'iO^Hm), (25) 

where F{z) = F{z) - z and F is defined by (HH]), i.e. F{z) = Kq + cA"1(z) - 
B{A~^{z)). Here Kq is a constant. In Fig. [1] (right) we plot the graph of a 
function F{z) for c > and /vq + c < 0. In such a situation, the function F has 
exactly two roots < z+ < 1/2 < z~ < 1 such that F{z^) = z^. Furthermore, 
F'{z~) > and F'(z+) < 0. As in the previous section [2.1.11 there exists, up 
to a shift in the argument ^, a unique solution z = z(^) connecting the steady 
states, i.e. lim^_>-|-oo z{£^) = z^. The corresponding traveling wave profile v = 
v{£,) given by v{£,) = A~^{z{^)) satisfies: 

lim w(^)=u±, v+ <{){£_)< i- , u'(^)<0, ^ £ M, 

where = A~^{z^), < v'^ < 1 < v~ . Again we can prescribe the limiting 
values < < 1 < and calculate the corresponding wave speed c > 
and constant Ko. Since = F{z^) = Ko + cv^ — B{v^) we have c = {B{v~ ) — 
B{v^))/{v^ — {;+) and A'o = —cv^ + B{v^). Since B = B + A we have 

-fl - v+]^ + 1 — 

S ^ L -■ (26) 

Summarizing the results of this section we conclude the following theorem. 

Theorem 1 For any limiting values < < I < there exists a speed 
c > given by j j7| ) such that the Hamilton-Jacobi-Bellman equation (0) with 
a range bound {9 < 1} has a solution V{x,t) such that the optimal response 
function 6*{x,t) given by (0) has the form of 0* — min(l, where ip{x,t) ~ 
v{x -\- c{T — t)) and viX) is a smooth and strictly decreasing function, 

lim^^ioo w(0 = v"^- 

The statement remains true if we consider the HJB equation hl8\) for 
the reduced optimal portfolio selection problem with the traveling wave pro- 
file (p{x,t) = v{x + c{T — t)), limj_^-j-oo '5(^) = and the optimal response 
function 0* = min(l, 1/0). Here the wave speed c> is given by h26\) for any 
prescribed limiting values Q < < 1 < . 

Proposition 1 For HJB equations (0) and M8\) there is no traveling wave 
profile with negative wave speed. 

Proof Indeed, for the traveling wave speed we have c = {B{v^ ) — B{v^)) / {v^ — 
v^) > because the function B is increasing. Analogously, c>QasB = B + A 
is an increasing function. Moreover, the profile v{^) is always a decreasing 
function, as there are no roots < z~ < ^ < 2;"'" < 1 of such that F'{z~) > 
and F'(z+) < 0. The latter follows from the inequality F(0) = A'o + f > 
Kq + c = F{^) and F{1^) = —00 in the case c < 0. The same statement holds 
true for the profile v. 
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Fig. 2 A graph of the function z (left) and the traveling wave profile v (right) for the 
parameter values = 1, d+ = 0.5, v~ = 2. 




Fig. 3 Graphs of the function 4){x,t) (left) and the response function 9*(x,t) (right). 

3 A HJB equation with more general drift and volatility functions 

One disadvantage of the HJB equations with a bounded constraint discussed 
in sections 12.11 and 12.21 is the fact that there is no travehng wave solution with 
increasing wave profile v. As a consequence, the optimal response function 
d*{x,t) is nondecreasing with respect to x. In this section we analyze a HJB 
equation with more general drift and volatility functions. Under suitable as- 
sumptions made on the model parameters, we shall prove that there exists a 
traveling wave solution having an increasing wave profile v and negative wave 
speed c < 0. 

We shall assume that a stochastic process follows a Brownian motion 
dXf = i/dt + a'^dWt 
with drift ji^ and volatility given by 

m 

where a,l3,uj,m G R,m > 1, and a; > arc model parameters. We again 
restrict our response strategy by 1 from above. The corresponding HJB 
equation ^ for the value function V{x,t) reads as follows: 



Traveling wave solutions to a Hamilton-Jacobi-BcUman equation 



9 



BV c 1 rP-V hV -1 

+sup|(«2 + + (/3 + c.0)^(x,i)} = 0, 

Ot e<i I. m ox ox J 

V{x,T) = u{x). (27) 

The simplified model discussed in the previous section corresponds to the 
choice a = /3 = 0, and m = 2. 

Again, supposing d^V > and d^V < the unconstrained optimal re- 
sponse strategy at (x, t) is the unique argument of the maximum of the 
function 

e ^ (a^ + -\9\"')dlV + + uj9)d^V. 
m 

Therefore 

where we have again employed the new variable defined by means of the 
Riccati-like transformation: 

uj dxV(x,t) 

Now, if (p{x,t) > 1 at {x,t) then d{x,t) < 1 and therefore, for the optimal 
response 6*{x,t), we have 9*{x,t) = 9{x,t). Hence the value function V{x,t) 
satisfies 

^ + (a- + -^---)^ + (/? + -v^---) - 0. 

Notice the following recurrent relations: 

dlV = -uj^d^V, dlV = (wV^ - ^dxH^)dxy- (28) 

Using the relation for d'^jV we can rewrite the equation for V in the form: 

dV dV ^ , ^ , rn-1 i_ 

— — 1-0-^=0, where g = p — uja ip + co lo "-i. (29) 

ot ox m 

Since 

dip ^ 1 dl,,v 1 dl,vdl,v ^ 1 dl,,v dl,v 

we obtain from and (|28p the equation 

dip Id^g 9 

which can be further rewritten in terms of the function ip as follows: 

dip , d^ f 2 m — 1 _ 1 \ ( Q 22,™^1 



dt dx^ V Ti J dx \ m / 

(32) 



10 



Naoyuki ISHIMURA, Daniel SEVCOVIC 



provided that ip{x,t) > 1. 

On the other hand, if (p{x, t) < 1 at (x, t) then 6{x, t) > 1 and therefore 
for the optimal response 9*{x,t) which is restricted by 1 from above we have 
d*{x,t) = 1. Then the value function V{x,t) satisfies the following equation: 

In view of relations (1251) we have 



dV ^dV , . ^ . 2 1 ^ 

— +5— =0, where g = fi + uj - uj(a H )tp. 

at ox m 

Using (PH)) . (|28l) and (PT|) (with g replaced by g) we finally obtain the following 
reaction diffusion equation 



dt 



«- + l)^ + ^((/3 + .).-.(.^ + lv^)=0, (33) 



which is fulfilled by if in the case when ip{x, t) < I. 

Similar to the simplified problem (a = /3 = 0, and m = 2) in both cases 
tp > 1 and tp < 1 we obtain that the function is a solution to (jlip . that is: 

where 

f (a' + ;7r)<^, for ^<1, 

^(^) = <^ (35) 

for ^>1, 

r (/? + C.)^ - Ujia^ + i^)^2 _ ^ni_l^ fo^ ^ < 1^ 

S(^) = <^ ^ ^ (36) 

[ - waV^ + ^'^('/'^ - 1), for > I- 

Notice that we have added appropriate constants to the definition of B in 
order to make it continuous across the value ip ^ I. Both functions A and B 
are continuous functions for (p > 0. Moreover, the function A is strictly 
increasing. 



3.1 A traveling wave solution with negative wave speed 

The aim of this subsection is to show, for suitable model parameters a,/3,uj 
and m, that there exists a traveling wave solution ip to (jlip with increasing 
wave profile v. To this end, we search for a traveling wave solution to (jlip of 
the form 

ip{x,t) = v{x + c{T -t)), a; e M, t e [0,r], 

with negative wave speed c < and where w = u(^) is a function defined on 
R. Inserting the traveling wave form of a solution p into (jlip we obtain that the 
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Fig. 4 A graph of the functions B{v) and Kq + cv for model parameters: m = 3/2, a = 
/3 = 0, c = -Ko = -0.1, tj = 1. 

function v = v{^) should fulfill identity (jl3p . Again, we introduce an auxiliary 
function z(^) = A{v{£,)). With this substitution, (p{x,t) = v{x + c{T — t)) is 
a traveling wave solution to pT|) if and only if the function z is a solution to 
the ODE: 

z'{0 = where F{z) = G{A~^{z)), G{v) = + cv - B{v). 

(37) 

Notice that the function G is G^ continuous for v > 0. The function A{v) is 
increasing for v > and its range is the interval (0, 1) if a = 0, and (0, oo) 
if a 7^ 0. In order to construct an increasing traveling wave profile v such 
that limt^ioo v{£,) = where < v~ < 1 < v~^, we first have to find roots 

= A{v^) of the function F such that F{z) > for < 2 < z+ and 
F'{z~) > 0,F'{z^) < 0. To this end, it is sufficient to investigate the roots 
and behavior of the function G{v) for positive v > 0. Clearly, if m > 2 then 
the function B is concave, B"{v) < 0,v 7^ 1. Thus G is convex and there are 
no roots such that G{v) > for v~ < v < . For this reason we cannot 
find the desired form of a traveling wave solution Lp for the case m > 2. 

The function v 1— > B{v) is strictly concave for v < 1. It has a local maximum 
at v~ = 2uj(a^i/m) ^'^'^ strictly convex for v > 1 provided that u > Q,a — 
and 1 < m < 2. Notice that, for a = /? = 0, 1 < m < 2 the maximum of 
B is attained at v_~ ~ m/2. In what follows, we restrict ourselves to the case 
a ~ (3 — 0, and m G (1, 2). Inspecting the behavior of the function B (see also 
Fig. m we conclude the following result. 

Theorem 2 Assume that a = /3 = 0, 7ti G (1,2), and oj > 0. Let the limits 
be such that v~ < < 1 and w+ > 1;+ > 1 where ~ m/2 < 1 and w+ > 1 
is the unique root of the secant equation B{v^) — B{v^) ~ B' {v^){v^ — v^). 

1. there exists a traveling wave speed c < given by 

B(v+) -B(v-) 

^ = + = 

u+ — v 

and intercept Kq = B{v~) — cv~ such that G{v^) = and G'{v~) > and 
G'{v+) < 0; 

2. there exists a solution V{x^t) to HJB with the optimal response func- 
tion 9*{x, t) given by 9* = min(l, ip^ ) where ip{x, t) ~ v{x + c{T — t)) 
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and v{£,) is a smooth and strictly increasing function, liin^_y±oo v{^) = 

Proof Part 1) follows from the behavior of the function B given by ([55]) for 
parameter values a = /? = 0, m G (1,2) and a; > 0. Notice that the condition 

> where t;+ is the unique root of the secant equation B{v^) — = 
B'{v~){v^ — v^), is necessary in order to find two roots of the equation Kq + 
cv — B{v) where c > B'{v~) with the property G{v~) = c — B'{v^) > and 
Giv+) = c-B'{v+) < 0. 

Part 2). Since the function A is strictly increasing and G'{v~) > and 
G'{v+) < the ODE ^ has a solution z = z{^) such that = and 

F{z) > for z- < z< z+ and -F'(z") > 0,F'{z+) < where z^ = 

Then there exists a traveling wave solution ip to of the form ip{x, t) = 
v{x + c{T — t)) with negative wave speed c < where v{£) = A~^{z{£)). 
Moreover, lim^^±oo "^(0 = = A~^{z^). If tp{x,t) > 1 then the optimal 
response 9*{x,t) is given by 6*{x,t) = 9{x,t) = ip{x, t)~ "-i . On the other 
hand, 6*{x,t) = 1 provided that (p{x,t) < 1, which concludes the proof. 

By a non-monotone traveling wave solution we mean a non-constant solu- 
tion whose derivative changes the sign several times. 

Proposition 2 There are no non-monotone traveling wave solutions to the 
fully nonlinear equations ill]) and ^34^ analyzed in sections \2.1.1\ and \3.1[ 
respectively. 

Proof We follow the ideas from [Q. If we assume to the contrary z'(Co) = 
then, as a consequence of the uniqueness of solutions to ODEs (|57|) and (jl4p . 
we obtain z{S) = const for all ^ G R. Consequently, the profile v is constant, 
a contradiction. 

We finish this section with two computational examples. In Fig. [5] we plot 
a graph of the function G{v) for the model parameters: m = 3/2, a = /3 = 
0, c = —Kq = —0.1 and a; = 1. In this case we have ii" = 1 is a root of G. 
Moreover, B'{1) = and so G"(l) = c - Hence G'{v-) > if and 

only if c > —u)^^—^. The function G also has a root v'^ = 10/3 > 1. Since 
TO = 3/2 and Kq + c = we have the analytic expression for w+ = > 1. 
Moreover, G'{v+) c + 3c2 < for w = 1, c -Kq = -0.1, and m = 3/2. 
In Fig. [S]we plot a graph of the function z and the traveling wave profile v. 
The functions (f>{x,t) and the optimal response function 0*{x,t) are depicted 
in Fig. [3 

For parameter values: to = 3/2, a = /3 = 0, c = —0.08, Kq = 0.1, and oj = 1 
the function has three roots. An unstable root is located at v~ « 0.888 < 1 and 
a stable root exists at « 4.488 > 1. Therefore the traveling wave solution 
ip(x,t) becomes strictly less than 1 for x — >■ — oo. Hence the optimal response 
function 9*{x,t) = min(l, (p~ "-i (.t, t)) is identically 1 for all sufficiently large 
negative x (see Fig. [5]). 

Finally, we give a justification of the assumption made on the value function 
V. In what follows, we shall prove that d^V < and dxV > 0. 
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Fig. 5 A graph of the function G{v) for model parameters: m = 3/2, a = /3 = 0, c = 
—Ko = —0.l,uj = 1 (left). In this case v~ = l,v~^ = 10/3. A graph of the function 
F{z) = G{A-^(z)) (right) with roots z' = 2/3, z+ ^ 0.97. 
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Fig. 6 A graph of the function z (left) and the traveling wave profile v (right) for the 
parameter values m = 3/2, a = /3 = 0,c = —Kq = — 0.1,aj = 1. 




Fig. 7 Graphs of <j){x,t) (left) and the optimal response function 8*{x,t) (right). 

Proposition 3 Suppose that the terminal condition V{x, T) = u{x) is a smooth 
function and there exist constants > such that X~ < —u"{x)/u'{x) < 
for all X G M. Then X~ /uj < ip(x, t) < /uj for all x £M. and t G [0, T] where 
if is a solution to equation 0j or il8\) and satisfying the terminal condition 
ip{x,T) = -{l/uj)u"{x)/u'{x), a; G M. 

Furthermore, for the value function V it holds: d^V{x, t) < and dxV{x, t) > 
for allxGR and t e [0,r] 

Proof Equations ([S]) as well (fTS)) are strictly parabolic partial differential equa- 
tions with a strictly positive diffusion coefficient A'((f>) which is uniformly 
bounded from below and above. Using the parabolic maximum principle we 
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Fig. 8 A graph of the optimal response function 8''{x,t) attaining the prescribed boundary 
9 = 1 (left), a graph of the optimal value 8*{x,t) computed by the stochastic dynamic 
optimization model [S] (right). 



conclude A /a; < (p{x,t) < A+/aj for all a; e M and t G [0,T], provided that 
X-/lu < (p{x,T) = -{1/lu)u"{x)/u'{x) < X+/uj for all a; e R. 

Concavity and the monotonicity of the value function V can be now de- 
duced from the fact -dlV{x,t)/dxV{x,t) = aj(/)(x,t) S (A~, A+). 



4 A motivation for analyzing HJB equations with range bounds 

The optimal allocation problem has a long history of research and much 
progress has been made so far (see for instance Dupacova [5]). For a com- 
prehensive overview of the stochastic dynamic optimization problems of the 
form ^ with the prescribed terminal value we refer the reader to papers by 
Merton [TT], Browne [7], Bodie et al. [5] and to the wide range of literature 
referenced therein. 

As an example of a process of the form ([T|) one can consider a stochastic 
process representing returns on the accumulated sum of a saver's portfolio 
consisting of volatile stocks and less volatile bonds. The time discrete version 
of the stochastic dynamic accumulation model has been proposed and analyzed 
by Kilianova et al. in [5]. In the time continuous limit, according to the pension 
savings accumulation model analyzed by Macova and Sevcovic in [9] , we have 
Xt = In Yt where the stochastic variable Yt represents the logarithm of a ratio 
Yt of the accumulated sum in the pension account and the yearly salary of a 
future pensioner. More precisely, by choosing a portfolio consisting of e [0, 1] 
part of stocks and 1 — 6 part of bonds one can derive the SDE for the ratio Yt : 

= (e + ^i'>Yt)t + a'^YtWt, (38) 

where / = efis + {l-0)tib and (cr*)^ = e^{as)^ + il-ef{ab)^+2ge{l-e)crsab 
(c.f. inj). In this model, ^J.s,^J'b and Os^Ob denote the expected returns and 
volatilities on stocks and bonds, respectively. Here g G ( — 1, 1) is a correlation 
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between returns on stocks and bonds. Using Ito's lemma for the variable Xt — 
InYt, we obtain the following SDE: 

Xf=(£e-^^*+/-i(a'^)^)t + a«W, (39) 

In a stylized financial market, one can assume the expected return fif, on 
bonds to be very small when compared to fis and also ab cts- Furthermore, 
we notice that the yearly contribution rate e > can be considered as a small 
model parameter (the value e = 0.09 was accepted in Slovakia, e = 0.03 was 
proposed in Czech republic, and, e = 0.14 was accepted in Bulgarian pension 
fund system). The asymptotic expansions of a solution to the corresponding 
HJB equation with respect to the small parameter < e <C 1 have been 
analyzed in [5]. By taking the formal asymptotic limit e ^ 0, 0, erf, —> 

and rescaling the time t n- t/a^, we end up with the following reduced SDE 
governing the variable Xt: 

x't^iujo-h^jt + ewt, (40) 

where uj = Hs/c^l is the Sharpe (reward-to- variability) ratio of returns on 
bonds. Such an underlying stochastic process corresponds to the one analyzed 
in Theorem [TJ The wave speed of the traveling wave solution from Theorem [1] 
is positive, c > and the profile ?; is a decreasing function, v' < 0. As a 
consequence, we have dtf > Q and dx^p < 0. Hence, for the optimal response 
function 9* = min(l, 1/ip), we obtain the following properties: 

dt9*{x,t)<0, dxO*{x,t)>Q, a; e M, i e (0,r). 

From the optimal allocation portfolio problem point of view, the decreasing 
behavior of the proportion of stocks 6* in the portfolio with respect to time 
can be interpreted as the necessity to be more conservative in investing when 
the time horizon T is approached. On the other hand, increasing dependence 
of 6* with respect to the accumulated property x indicates an increase of 
the investor's risk selection preference with the increase of the accumulated 
property x. 

5 Conclusions 

In this paper we have investigated traveling wave solutions to a fully nonlinear 
parabolic PDE arising from the model of a stochastic dynamic optimal allo- 
cation problem. Our approach is based on the method solving and analyzing 
the fully nonlinear parabolic partial differential equation which can be con- 
structed from the corresponding Hamilton-Jacobi-Bellman equation. Under 
certain assumptions made on the form of the underlying stochastic process, 
we found a range of parameters for which a traveling wave solution can be 
constructed. More precisely, we identified parameters for which the traveling 
wave profile is monotonically decreasing or increasing and it has positive or 
negative wave speed, respectively. 
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